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Ab-initio density functional theory calculations are performed to study the electronic properties 
of a M0S2 monolayer deposited over a SiCb substrate in the presence of interface impurities and 
defects. When M0S2 is placed on a defect-free substrate the oxide plays an insignificant role, since 
the conduction band top and the valence band minimum of M0S2 are located approximately in the 
middle of the SiCb band-gap. However, if Na impurities and O dangling bonds are introduced at the 
SiC>2 surface, these lead to localized states, which modulate the conductivity of the M0S2 monolayer 
from n- to p-type. Our results show that the conductive properties of M0S2 deposited on SiCb are 
mainly determined by the detailed structure of the MoS2/Si02 interface, and suggest that doping 
the substrate can represent a viable strategy for engineering MoS2-based devices. 

PACS numbers: 



I. INTRODUCTION 

Recently MoS2-based layered transition metal chalco- 
genides (LTMDs) have attracted considerable attention 
due to their potential for constructing low dimensional 
nano-structures for a variety of application^^. The elec- 
tronic properties of M0S2 show a strong dependence on 
thickness^, i.e. on the number of atomic layers forming 
a given sample. In particular, M0S2 monolayers, which 
display a substantial direct band-gap, represent a semi- 
conducting alternative to graphene, which is a metal in 
its pristine form. Although the band-gap of graphene can 
be opened by fabricatingnanoribbons^ or by depositing 
it on a suitable substrate 6 , this comes to the prize of de- 
teriorating, in a somehow uncontrollable way, the carrier 
mobility due to edges and impurity scattering 7 . In con- 
trast, the low dimensionality, the small amount of dan- 
gling bonds, and their typical high crystalline form, make 
the performances of LTMD-based transistors comparable 
to those of existing Si-based ones^H^. In particular tran- 
sistors made from M0S2 monolayers have been recently 
fabricated, showing a mobility of at least 200 cm 2 /V-s at 
room temperature, an on/off current ratio of 10 8 and low 
standby power dissipation^. 

Interestingly, both n-type^M^ and p-type^LL^ conduc- 
tivities have been reported in ultra-thin M0S2 layers de- 
posited on Si02- The conducting behavior of M0S2 there- 
fore seems to depend on the experimental details and 
an explanation for the specific current polarity (n- or 
p-type) remains far from being clear. Note that no in- 
tentional doping was introduced in the experiments men- 
tioned above, so that the source of the different carrier 
types should be intrinsic to the M0S2 layer, to the sub- 
strate and to the interaction between the two. 

The possible creation of Mo and/or S vacancies dur- 
ing the growth cannot be the cause of the various con- 
ductive properties, since vacancies create deep levels at 
mid-gap in the bandstructure of M0S2 monolayers 17 . No- 
tably, disorder at the semiconductor /substrate interface 
in general plays a crucial role in determining the con- 
ductive properties of ultra- thin devices. For example, 



for GaAs nanowires it has been demonstrated that upon 
decreasing the nanowire diameter the interface-mediated 
conductivity gradually becomes dominant over the bulk 
one^. Such surface sensitivity can also be used to one's 
advantage. For instance an ambipolar transistor has been 
realized in M0S2 thin flakes by contact with a ionic liquid 
environment!^, which as well affects the interface prop- 
erties. Since M0S2 monolayers are placed on insulators 
in practically any device architecture, it is important to 
identify the possible effects that the substrate has on the 
conductivity. 

The defects responsible for the conductive properties 
of low dimensional devices are expected to be extrin- 
sic in nature, such as charged impurities at the inter- 
face between the conductive channel and the substrate. 
These lead to an inhomogeneous Coulomb potential for 
both conduction and valence band electrons. Such charge 
traps have been identified to be in the form of adsorbates 
or defects at the su rface of the underlying substrate in 
the case of graphen d 20 * 21 *. Likewise, temperature depen- 
dent transport measurements on thin M0S2 layers, down 
to the monolayer limit, suggest that trapped charges at 
the Si02 surface could be responsible for the observed 
n-type behavior, when M0S2 is deposited on SiOsr^. 

In general when charged traps are located at an inter- 
face, they influence the depletion/accumulation of elec- 
trons in the conducting channel up to a certain thickness, 
which is proportional to the channel screening length. 
This distance depends on different physical features, such 
as the nature and the density of the traps, and the 
electronic properties of the channel. For conventional 
semiconductors it typically reaches up to a few nanome- 
ters. For instance it has been recently demonstrated that 
charged trap states at the substrate/channel interface sig- 
nificantly affect the conductivity of GaAs nanowires up 
to diameters of about 40-70 nm 18 . More dramatic effects 
are expected for layered compounds down to the single 
layer limit, in which essentially all the atoms are at the 
interface with the substrate and the channel vertical di- 
mension is certainly shorter than the screening length. 
Recently a reduction in conductivity with increasing the 
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M0S2 film thickness has been observed in MoS2-based 
transistors, where Si02 was used as back gate^. This 
suggests that for the MoS2/Si02 system the transport is 
interface-mediated, as intrinsic defects, homogeneously 
distributed in M0S2, would not lead to any dependence 
of the conductivity on thickness. 

In order to shed some light on the effects that trap 
states at the Si02 surface have on the conductive prop- 
erties of MoS2/Si02 hybrid systems, we have performed 
state of the art first principle electronic structure calcu- 
lations. In particular we have considered the case when 
the traps are due to impurities such as immobile Na and 
H atoms, and O-dangling bonds. The paper is organized 
as follows. In the next section we briefly describe our 
computational techniques and we provide details of the 
simulations performed. Then we proceed with presenting 
the results of this work in the context of recent experi- 
ments, and finally we conclude. 



II. METHODOLOGY 

In order to investigate the influence of a Si02 substrate 
on the electronic properties of a M0S2 monolayer, ab- 
initio calculations are pe rformed by using density func- 
tional theory (DFT/^El w ithin the generalized gradient 
approximation (GGA) of the exchange and correlation 
(XC) functional as introduced by Perdew, Burke and 
Ernzerhof (PBE) ^ and numerically implemented in the 
SIESTA codeP. In our calculations, a double-^ polar- 
ized^ numerical atomic orbital basis set is used for all 
the atoms and the Troullier-Martins scheme is employed 
for constructing norm-conserving pseudopotentials 29 4 An 
equivalent plane wave cutoff of 350 Ry is chosen in all the 
simulations and the Brillouin zone is sampled by using 
an equivalent /c-grid cutoff of 17 A. Relaxed geometries 
are obtained with the conjugate gradient method, where 
all the atoms in the supercell are allowed to relax until 
the force on each atom is less than 0.02 eV/A. 

A trap state is usually formed when an energy level as- 
sociated either to a defect or an impurity appears within 
the energy gap of the host material. Such trap states 
influence the charge transport properties mainly in two 
ways. Firstly, if the traps are charged, they will capture 
a hole or an electron from the environment. This pro- 
duces a modification of the electrostatic potential, which 
in turns shifts the level alignments in the system, and 
thus affects the conductivity. Secondly, they can also in- 
crease the carrier concentration and provide pathways for 
electrons or holes to hop. The efficiency of this process 
depends on the amount of localization of the states asso- 
ciated to the defect site. If the energy of the localized gap 
state is close to either the valence band maximum (VBM) 
or the conduction band minimum (CBM), then at a given 
temperature some of these charges will be transferred ei- 
ther to the conduction or to the valence band, where they 
may contribute to increase the system conductivity. 

Whether or not one can describe with ab initio calcu- 



lations such mechanisms depends crucially on the ability 
of computing accurately the energy levels of the system. 
The use of the GGA (or of the local density approxima- 
tion - LDA) for electronic structure calculations of defect 
levels is, in general, problematic. One reason is the typ- 
ical underestimation of the energy gap and the related 
incorrect alignment of the energy levels of hybrid sys- 
tems. For instance an artificially reduced band-gap may 
erroneously bring deep traps in reson ance with either the 
conduction or the valence bancP^U. A second source of 
error is the incorrect description of the charge localization 
at the defect site, a feature that usually leads to predict 
defects to be too shallow 3 -^. Atomic self-interaction cor- 
rection 

(ASICPH has been proved to overcome these 
deflciencie d 35 * 36 *. Therefore we also perform additional 
LDA+ASIC calculations to verify the robustness of the 
LDA/GGA results. In particular we set the ASIC scaling 
parameter to a = 0.5, a value which is generally appro- 
priate to mid-gap insulators 3 - 3 -^ 



III. RESULTS AND DISCUSSION 

A. Defect-free Si02 interface 

Substantial experimental efforts have been devoted to 
deposit ultra-thin M0S2 layers onto Si02 in order to 
demonstr ate tran sistor operation, down to the single 
layer limitP^^lii] Usually amorphous oxides are used as 
substrates. However, in order to avoid the computational 
complexity of a highly disordered structure, a crystalline 
Si02 substrate is simulated here instead. This also al- 
lows us to systematically determine the effects of indi- 
vidual defects and impurities on the electronic structure 
of a M0S2 layer. Our unit cell is constructed as a slab 
containing at least 6 Si atomic layers of a-quartz and an 
adsorbed M0S2 monolayer. At least 15 A of vacuum are 
included at the slab boundaries to avoid the spurious in- 
teraction between the slab periodic images. We consider 
the modified oxygen-terminated (0001) Si02 surface in 
order to simulate the most experimentally relevant con- 
ditions. 

Two primary structures for the oxygen-terminated 
Si02 (0001) surface are possible, depending on whether 
the termination is with the siloxane group (Si-O-Si) or 
with the silanol one (Si-OH). Both surfaces can form de- 
pending on the surface treatment^. The siloxane recon- 
struction at room temperature forms an O-terminated 
surface with an outermost six-membered ring structure, 
as shown in Fig.[TJa, b). Under annealing in ambient con- 
ditions it becomes hydroxylated (Si-OH) and the recon- 
struction transforms into the silanol one, which presents 
on the surface a zigzag H-bonded network [see Fig. [TJc, 
d)]. In both cases in our simulations the dangling bonds 
on the Si-terminated bottom surface are saturated by hy- 
drogen. 

The optimized lattice constants of the pristine Si02 
and M0S2 are 4.91 A and 3.19 A, respectively. We there- 
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FIG. 1: (Color online) Side and top views of reconstructed 
structures for the O-terminated SiCb (0001) surface (a and 
b), and the fully hydroxylated SiCb (0001) one (c and d) 
(color code: cyan — > Si, red — > O, violet — > H). 




FIG. 2: (Color online) Top view of the optimized structure 
of M0S2 placed on a defect-free (a) siloxane and (b) silanol 
surface (color code: light grey — > Mo, yellow — > S, cyan — > Si, 
red — >> O, violet — >> H). 



fore construct a hexagonal supercell in the plane, with a 
9.69 A-long side, so that the lattice mismatch between 
Si02 and M0S2 is minimized to ~1.2 %. The GGA cal- 
culated band-gap of Si02 and of a M0S2 monolayer are 
6.20 eV and 1.49 eV, respectively. The small strain ap- 
plied to the M0S2 monolayer changes only little the elec- 
tronic structure. The band-gap remains direct at the K 
point and it is only reduced by 0.22 eV from the value of 
1.71 eV obtained for the unstrained case. Similarly to the 
case of graphene 3 -^, we expect that the electronic struc- 
ture of a M0S2 monolayer is only marginally affected by 
its local arrangement on the Si02 substrate. Therefore, 
as a representative configuration we use the arrangement 
of Fig. [2j where an oxygen atom is situated at the hollow 
site of the Mo surface triangles. 

We start our discussion by presenting the properties 
of the defect-free hybrid MoS2/Si02 system. The equi- 
librium distances, do, between the Si02 and the M0S2 
surfaces are 3.01 A and 2.98 A for siloxane and silanol, re- 
spectively. Here do is denned as the vertical separa- 
tion between the top-most O layer in the Si02 surface 
and its nearest S layer in M0S2. These values are sim- 



ilar to the distance between two M0S2 monolayers that 
we calculate to be 3.17 A. The binding energy of the 
MoS 2 /Si0 2 system is given by E h = E MoS2 + ^si0 2 - 
^MoS 2 /Si0 2 > where E Mo s 2 , ^Si0 2 , and ^MoS 2 /Si0 2 are to- 
tal energies of the isolated M0S2, the isolated Si02 slab, 
and of the MoS2/Si02 hybrid system, respectively. We 
find E\> for siloxane and silanol to be respectively 0.14 eV 
and 0.16 eV per primitive M0S2 unit cell. These binding 
energies are close to that between two M0S2 layers (0.20 
eV/unit cell), which are bound together by the rather 
weak van der Waals forces. As such, our results show 
that M0S2 is weakly bound also to the Si02 surface, 
in agreement with recent experimental results that have 
measured the interaction between M0S2 and an underly- 
ing Si02 substrate to be negligible 3 -^. 

Note that in general GGA- type XC-functionals do not 
describe accurately van der Waals forces. However, it 
has been shown that the LDA/GGA is able to repro- 
duce the interlayer spacing and the binding energy of 
layered chalcogenides^. We have then verified that our 
calculated do for bulk M0S2, do = 3.08 A, is i n g ood 
agreement with the experimental value of 2.96 and 
also with the previously calculated theoretical estimate 
of 3.05 A^. Moreover, in order to take into account pos- 
sible small deviations of the relaxed distance from the 
experimental value due to the XC functional used, we 
have evaluated the electronic structure for do within the 
range do ±0.5 A, and we have found that the results 
change little with varying the distance. 



B. SiC>2/MoS2 composite with siloxane 
reconstruction 

We now move to study the electronic structure of a 
M0S2 monolayer deposited on Si02 by starting with the 
siloxane surface. In particular we consider first the situa- 
tion where Si02 is defect-free. Fig.[3Fc) shows the density 
of states (DOS) of the hybrid Si02/ M0S2 system, which 
remains semiconducting with a band-gap of 1.48 eV, i.e. 
with the same band-gap of a free-standing M0S2 mono- 
layer having the same lattice parameters. Both the va- 
lence and the conduction bands of the hybrid compound 
are associated to M0S2. We note that the projected DOS 
(PDOS) over M0S2 extends into the Si02 band-gap, and 
the total DOS of the combined material is essentially 
given by the superposition of the DOSs of the pristine 
slab of Si02 [Fig. [3^a)] and of the M0S2 monolayer 
[Fig. |3jb)]. Both the conduction and the valence bands 
of Si02 are located at least 1.5 eV away from those of 
M0S2. As a consequence, no charge transfer between the 
substrate and M0S2 occurs. Importantly, one of the ba- 
sic criteria for the selection of the gate oxide is fulfilled 
here, namely that the oxide should have a bands offset 
of over 1 eV for both the conduction and valence band 
in order to create a large barrier for both electrons and 
holes 43 . Our results show that the conductivity of M0S2 
is not influenced by the underlying defect-free Si02 sub- 
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FIG. 3: (Color online) Electronic structure of the Si02/MoS2 
hybrid system when various defects are present at the SiCb 
siloxane surface, (a) The total DOS for the defect-free sur- 
face (black, dashed curves), and when either Na (green, solid 
curves), or H adsorbed (magenta, solid curves) are adsorbed, 
(b) The DOS of a pristine free-standing M0S2 monolayer. The 
total DOS and the PDOS for MoS 2 , when the MoS 2 mono- 
layer is placed on the defect-free siloxane surface (c), on a 
siloxane surface with one adsorbed Na, or (d) with one ad- 
sorbed H. The blue dashed vertical line indicates the Fermi 
level, which has been set to zero in all the panels. The red 
shaded areas indicate the M0S2 PDOS. Positive and negative 
DOS are respectively for spin up (majority spins) and spin 
down (minority spins) electrons. 



strate. Therefore the measured n-type or p-type conduct- 
ing properties of M0S2 on Si02 must be due to defects 
and impurities. 

Localized states, arising from impurities or defects 
within the oxide substrate or at the interface with the 
conducting channel, can redefine the effective Fermi level 
of the hybrid system, as illustrated schematically in 
Fig.|4j Depending on the alignment of the gap states with 
respect to the M0S2 valence and conduction bands, the 
system can switch from n-type [see Fig. |4jb)] to p-type 
[see Fig. Qc)]. Therefore, such trap states are expected 
to give significant contributions to the conductivity of 
these low dimensional systems. In the layered structure 
considered in this work, the trap states are expected to 
be located at the interface between the LTMDs and the 
substrate, not in the LTMDs themselves, which usually 
are highly defect-free. Trap states at the Si02 surface 
can have a wide range of origins, such as immobile ionic 
charges, Si02 surface dangling bonds, and foreign impu- 
rities adsorbed on the surface 4 ^. In literature densities of 
trap states on Si02 are reported in the range 45 10 10 -10 14 
cm -2 . As representative dopants, here we consider two 
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FIG. 4: (Color online) Schematic band diagram for M0S2 
placed on the defect Si02 substrate (a), and on a substrate 
including a defect-induced donor (b) or acceptor (c) level. 
This demonstrates the modulation of the conductivity from 
n-type to p-type as the impurity state redefines the Fermi en- 
ergy in the oxide. The energy levels Ey , Ec, Ef and E vac 
define the valence band maximum (VBM), the conduction 
band minimum (CBM), the Fermi energy, and the vacuum 
level, respectively. The subscripts 1 and 2 refer to Si02 and 
M0S2, respectively. The blue dashed- lines indicate the Fermi 
energy of the hybrid system (common to Si02 and M0S2). 
The thick green line in (b) indicates the donor level and the 
thick red line in (c) represents the acceptor state in the oxide. 
Note that in general due to charge transfer from M0S2 to the 
gap states, and the related dipole formation, the level align- 
ment between Eyi and Ev2 will also change in the defective 
systems. 




FIG. 5: (Color online) The optimized geometry for M0S2 
placed on (a) the siloxane surface incorporating a Na impu- 
rity and (b) the dangling oxygen bond on the silanol surface, 
obtained by removing a H atom. The arrows indicate the po- 
sitions of the defects on the surface (Color code: green — > Na, 
while the other colors are the same atoms as in Fig. [2|. The 
arrows indicate the location of the impurities/defects. 



possible candidates: Na atoms and Si02 surface oxygen 
dangling bonds. 

During the synthesis and the sample preparation, Si02 
can adsorb relatively light impurities such as Na and 
at its surface. In order to simulate the effects of such im- 
purities on the electronic structure of the M0S2 channel 
a single Na atom is placed on top of the siloxane Si02 
surface. Given the lateral dimension of our supercell, 
this corresponds to an impurity density of ~10 14 cm -2 , 
which is close enough to the recently reported values of 
trap states densities, reaching up to ~10 13 cm -2 for thin 
M0S2 layers deposited on SiOj^^. The most energeti- 
cally favorable binding position for Na is found to be at 
the center of the surface oxygen triangle [see panel (a) of 
Fig. [5] . A Na adatom adsorbed on a pristine Si02 sur- 



face creates a deep donor state in the DOS, with a single 
particle level at about 2 eV below the Si02 CBM [see 
Fig. [3^a)]. Note that such state is singly occupied and 
therefore spin-splits in our spin-polarized calculations, 
with the empty minority spin state (spin down) laying 
approximately 1 eV below the CBM and 1 eV below the 
Fermi level. 

When a M0S2 monolayer is deposited over the the 
Na-doped Si02 surface, do increases to 3.24 A at the 
edges of our unit cell, whereas at the Na site the O-S 
distance becomes 3.45 A. The enlargement of the bind- 
ing distance compared to that of the pristine Si02/MoS2 
system is a direct consequence of the Na intercalation at 
the interface. The electronic structure of the compos- 
ite is strongly affected by the presence of the Na ion, 
as shown in Fig. [3^c). Also in this case the total DOS 
appears as a direct superposition of those of Si02 and 
M0S2. However the presence of the Na filled state shifts 
the Fermi level, which now gets pinned just below the 
MoS 2 CBM. The resulting DOS around E F is thus that 
of the defect-free M0S2 conduction band with the addi- 
tion of a Na-derived impurity level positioned below it. 
Hence, the gap state is moved below the Fermi energy, 
resulting in a very small activation energy for the transfer 
of electrons from Na to the M0S2 conduction band. This 
is the situation schematically presented in figure [4^b), 
which leads to n-doping. 

If we now replace Na with H on the Si02 siloxane sur- 
face, the associated filled gap state lies deep in the Si02 
band-gap [see Fig. [3|a)], despite the fact that H and Na 
share the same s-like valence. The same situation persists 
in the composite [see Fig.|3|e)], where the H-derived filled 
spin- up level remains at mid-gap, approximately 0.5 eV 
above the VBM, while the empty spin-down one is nearly 
resonant within the conduction band. This situation 
however does now lead to doping so that H can not influ- 
ence the conductivity of the MoS2/Si02 structure. The 
quantitative difference found between the results for the 
Na and the H case show that, in order to obtain n-type 
character, only impurities with rather small ionization 
potential are relevant. These can transfer one electron to 
the M0S2 conduction band with small activation energy. 
Such activation energy is a key factor in the determina- 
tion of the threshold voltage, Vth, required to operate a 
transistor in the on-state. As a consequence, the exper- 
imentally measured values for Vth, which show a large 
variation for different samples^, are then attributable to 
varying concentrations and properties of the trap states 
from sample to sample. 



C. Si02/MoS2 composite with silanol 
reconstruction 



Next we move to examine the case of the Si02 sur- 
face with silanol reconstruction, whose DOS is presented 
in Fig. [6^a). Similarly to the siloxane case, the PDOS 
for the defect-free MoS2/Si02 composite [see Fig. [ljb)] 
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FIG. 6: (Color online) Electronic structure of the Si02/MoS2 
hybrid system when various defects are present at the Si02 
silanol-reconstructed surface, (a) The DOS for the defect- 
free surface (black, dashed curve) and for the one where one 
O dangling bond is induced by a single H removal (green, solid 
curve labelled as SiO*). The total DOS and the M0S2 PDOS 
for the Si02/MoS2 composite when the M0S2 monolayer is 
placed on (b) the defect-free surface, (c) on the surface with 
a single adsorbed Na atom, and (d) on the surface with a 
O dangling bond created by removing a single H atom. The 
blue dashed line indicates the Fermi energy, which is set to 
zero in all panels. The red shaded areas indicate the M0S2 
PDOS. Positive and negative DOS are respectively for spin 
up (majority spins) and spin down (minority spins) electrons. 



corresponds to a superposition of the DOSs of the iso- 
lated MoS 2 [Fig. [3fb)] and Si0 2 [Fig.^a)] components, 
indicating weak interaction between the two materials. 
When a Na atom is intercalated between the silanol sur- 
face and the M0S2 layer, we find that the system becomes 
n-type [Fig.^c)], in the same way as for the siloxane sur- 
face. This indicates that Na is an efficient n-dopant for 
M0S2 on Si02 regardless of the surface reconstruction. 

In general thermal annealing of the silanol surface cre- 
ates under-coordinated oxygen atoms (Si-O*). These ap- 
pear as stable surface defect centers and act as typical 
charge traps in oxygen rich SiOj^, since they are able 
to capture an extra electron in their dangling bond. In 
our calculations, such defects are created on the Si-OH 
surface by removing a H atom from the top surface [see 
Fig. [5jb)] . For such a defect we find that the empty 
acceptor state is created ~0.9 eV above the Si02 VBM 
[see Fig. [6ja)]. Once M0S2 is layered onto the surface, 
the value of do at the boundary of our H-deficient unit 
cell is do — 2.98 A, which is approximately equal to that 
for the pristine surface, whereas at the dangling bond 
site the O-S distance is significantly reduced to 2.68 A. 
When placing the M0S2 monolayer on this defective sur- 
face, the dangling bond state gets filled by capturing an 
electron from the M0S2 valence band, so that the Fermi 
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FIG. 7: (Color online) Density of states for the defective 
Si02/MoS2 composite calculated with the ASIC XC func- 
tional. In panel (a) we report the DOS for the siloxane re- 
construction with an intercalated Na atom [corresponding to 
Fig|3jd)], while in (b) that for the silanol reconstruction and 
an O dangling bond obtained by removing a surface H atom 
[corresponding to|6jd)]. 

energy now lies just below the M0S2 VBM [see Fig.^d)]. 
This is the level alignment presented in Fig. [4jc), which 
makes the composite p-type. Note that the rather high 
density of oxygen dangling bonds in our system causes a 
large surface charge density dipole, which shifts the M0S2 
DOS downwards in energy by more than 1 eV with re- 
spect to the Si02 substrate. By modulating the density 
of such defect types one may be able to change such a 
shift. 



D. Robustness of the results against the choice of 
XC functional: ASIC 



portantly the ASIC result is rather close to that calcu- 
lated^ with the hybrid Heyd-Scuseria-Ernzerhof (HSE) 
exchange-correlation functional. This is however larger 
than the optical band-gap of 1.90 eV measured experi- 
mentally for M0S2 monolayers^. The apparent contra- 
diction can be solved by noting that the optical exci- 
tations involve excitons with a large binding energy of 
the order of 1 eV, as confirmed by many-body calcula- 
tions^. Thus, one expects that the true quasi-particle 
spectrum has a band-gap of approximately 1.9+1=2.9 
eV, in good agreement with that computed with the 
GW scheme, either at the first order level 49 (2.82 eV) 
or self-consistently 51 (2.76 eV). As such, the ASIC de- 
scribes M0S2 with a band-gap larger than the measured 
optical one and it provides an improved description over 
that of the GGA. 

We now go back to the Si02/MoS2 composite and in 
Fig. [7] we report two representative results for the case 
of Na adsorbed on the siloxane surface and for that of 
the O dangling bond on the silanol one. We find that 
for the first case, although the band-gaps of the two 
parental materials are both increased, the Fermi energy is 
still pinned at the bottom of the M0S2 conduction band 
[Fig. (7ja)]. As a consequence Na still leads to a n-type 
semiconducting character with small activation barrier. 
Similarly the O dangling bond on the silanol-terminated 
surface leads to a p-type semiconducting character [see 
Fig. (7^b)], with the Fermi energy positioned below the 
M0S2 valence band. This indicates that our two main re- 
sults remain unchanged whether calculated at the GGA 
or ASIC level, i.e. they are robust with respect to the 
choice of exchange-correlation functional. 



IV. CONCLUSION 



Finally, in order to verify that the calculated level 
alignment is robust against the choice of exchange and 
correlation functional, we have repeated our calculations 
by using the ASIC scheme. As expected the ASIC func- 
tional increases the band-gap of M0S2 and Si02 respec- 
tively to 1.73 eV and 8.02 eV (for the same strained hy- 
brid structure used in the previous sections). In the case 
of Si02 this brings the calculated value sensibly closer to 
the experimental one of 8.9 eV 48 , as expected from the 
ASIC when dealing with an insulator whose va lence and 
conduction bands have different orbital content ! 33 ^ 34 !. 

The situation for M0S2 is more complicated and de- 
serves a detailed discussion. In this case the band-gap is 
defined by bands dominated mainly by Mo-d orbitals and 
the ASIC opens it only marginally. For a free-standing 
M0S2 monolayer the ASIC (a = 0.5) returns a direct 
band-gap of 2.03 eV (compared to a GGA gap of 1.71 eV). 
Note that the LDA value is 1.87 eV so that the LDA al- 
ready partially opens the gap with respect to the GGA. 
It is also notable that an enhancement of the screening 
parameter a to a = 1 (full atomic correction) produces 
a marginal further increase of the gap to 2.10 eV. Im- 



The effects of the Si02 substrate on the conductiv- 
ity of a semiconducting M0S2 monolayer are investigated 
with first principles density functional theory calcula- 
tions. The defect-free Si02 surface does not affect sig- 
nificantly the electronic properties of M0S2 due to their 
weak mutual interaction. As such the conductive proper- 
ties of M0S2 do not change and Si02 appear as an ideal 
gate material. However, when Na atoms are placed at 
the Si02/MoS2 interface, a shallow donor trap state is 
created just below the CBM of the hybrid Si0 2 /MoS 2 
composite. The small activation energy makes the hy- 
brid MoS2/Na-Si02 system a n-type semiconductor even 
for rather low temperatures. Interestingly, the behavior 
is different for H adsorption, where the impurity level is 
created ^0.9 eV below the CBM, resulting in a stable 
localized charge that cannot be easily promoted to the 
CBM and does therefore not affect the conductivity. 

In contrast, in the case of oxygen dangling bonds on the 
silanol-terminated Si02 surface, the Fermi energy of the 
MoS2/Si02 system is located just below the VBM, mak- 
ing the system a p-type semiconductor. These results 
shows that the conductivity of ultra-thin semiconduct- 
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ing LTMDs changes from n-type to p-type depending on 
the charge-polarity of the traps, as well the energy level 
alignment of the trap states within the LTMDs band gap. 
These kind of trap states at the SiC>2 surface are likely to 
be at the origin of the observed change in conductance in 
different experimentally realized MoS2-based transistors. 
Intriguingly, our results suggest the possibility of inten- 
tionally doping M0S2 by depositing different adsorbates 
over the substrate Si02 surface. This can pave the way 
for a new strategy in the design of two-dimensional de- 
vices, where the electronic properties of the channel are 
engineered by manipulating those of the substrate. 
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